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INTRODUCTION 

Since extinction has been observed in an oscillatory manner in Le> 1 premixed flames [1], it is 
not unreasonable to expect that extinction could occur in an unsteady manner for diffusion 
flames. Indeed, near-limit oscillations have been observed experimentally under microgravity 
conditions for both candle flames [2] and droplet flames [3], Furthermore, the analysis of 
Cheatham and Matalon [4] on the unsteady behavior of diffusion flames with heat loss, identified 
an oscillatory regime which could be triggered by either a sufficiently large Lewis number (even 
without heat loss) or an appreciable heat loss (even for Le= 1). 

In light of these recent understanding, the present investigation aims to provide a well- 
controlled experiment that can unambiguously demonstrate the oscillation of diffusion flames 
near both the transport- and radiation-induced limits. That is, since candle and jet flames are 
stabilized through flame segments that are fundamentally premixed in nature, and since premixed 
flames are prone to oscillate, there is the possibility that the observed oscillation of these bulk 
diffusion flames could be triggered and sustained by the oscillation of the premixed flame 
segments. Concerning the observed oscillatory droplet extinction, it is well-known that gas-phase 
oscillation in heterogeneous burning can be induced by and is thereby coupled with condensed- 
phase unsteadiness. Consequently, a convincing experiment on diffusion flame oscillation must 
exclude any ingredients of premixed flames and other sources that may either oscillate 
themselves or promote the oscillation of the diffusion flame. The present experiment on burner- 
generated spherical flames with a constant reactant supply endeavored to accomplish this goal. 
The results are further compared with those from computational simulation for further 
understanding and quantification of the flame dynamics and extinction. 

SPECIFICATION OF EXPERIMENT AND COMPUTATION 

An inverse flame configuration was employed with the oxidizer being ejected from a porous 
sphere burner ( 20pm pore-diameter, 1.27cm sphere diameter) into a low-density fuel 
environment comprising of He and H 2 at a pressure of 0.096atm to effectively minimize 
buoyancy. Various inert gasses (N 2 , CO 2 , He) were used to dilute the oxidizer in order to change 
the Lewis number and radiative properties of the mixture so that both the transport-induced and 
radiative-induced limit could be achieved. A Xybion IMC-201 intensified multi- spectral video 
camera with a UV-transmissive lens was adopted for visualization. The excited CO 2 (CO 2 *) 
chemiluminescence was used as a marker for detecting the CCU-diluted flames, whereas OH* 
due to transition L — » HI, observable at 305. 4nm, was used for the N 2 -diluted flames. The 
radiometer was placed 14 cm from the center of the burner to detect emission changes in the 
flame before extinction. 

Extinction was triggered by gradually decreasing the H 2 mole fraction in the ambient. Four 
oxidizer mixtures, hereafter referred to as cases (a), (b), (c) and (d), diluted with different 
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combinations of CO 2 , No, and/or He flow rates were selected to modify the transport and 
radiative properties of the mixture. Table 1 lists their respective mixture composition and mass 
flow rates. 

The spherically symmetric unsteady diffusion flame was simulated using a modified version 
of the code developed by Ref. [5] where terms representing unsteadiness, radiation, heat loss, and 
boundary conditions for non-premixed burning were added. The flame was assumed to be 
optically thin with COo, CO, and H 2 O as radiating species, with respective emission/absorption 
coefficients obtained from Ju et al. [6]. The mass flow rate m, the burner temperature T s (432K), 
the oxidizer and ambient fuel compositions, and the chamber pressure p were selected to 
simulate the experiment. The size of the computational domain was selected to be the effective 
radial distance (24.8cm) from the burner surface to the chamber wall. The boundary condition at 
the end of the computational domain on temperature was the chamber wall temperature measured 
in the experiment (7)^30510. The H 2 /O 2 reaction mechanism used was that of Mueller et al. [7]. 
To simulate the effect of adding CO 2 as diluent, Mueller’s mechanism was augmented with 9 
additional elementary reactions, describing the oxidation of CO and the formation/consumption 
of HCO from the moist CO mechanism of Kim et. al. [8]. Helium was also added to the 
mechanism simply by assuming the same third body efficiencies as argon. 

EXPERIMENTAL AND COMPUTATIONAL RESULTS 

In general, as the ambient H 2 mole fraction was decreased, the flame luminosity became 
weaker and the flame size increased. Furthermore, the occurrence of flame oscillation in terms of 
luminosity and flame size was not visually discernable. The experimentally determined ambient 
Ho mole fractions at flame extinction are listed in Table 1. Although all cases appeared to 
extinguish without pulsating visibly in terms of the flame size, the recorded emission data for 
cases (a) and (b) revealed that the flames pulsate at 2.1 and 2.3Hz, shown in Fig. 1, prior to 
extinction, as was predicted in Ref [9]. For cases (c) and (d), no oscillation was observed as 
shown in Fig. 2. 

The steady-state response of spherically symmetric diffusion flames was computationally 
simulated for the experimental cases (a) through (d). In general the maximum flame temperature 
decreased while its location increased as the ambient fuel concentration was decreased. When the 
maximum flame temperature and the fuel molar concentration were plotted, a typical extinction 
flame response curve was obtained. Whether the observed extinction limit is transport-induced, 
for which T max increases with increasing m, or radiation-induced, for which T max decreases with 
increasing m, may be determined by perturbing the mass flow rate from the burner, as discussed 
in Ref. [9]. For the four mixtures investigated here, this interrogation reveals that cases (a) and 
(b) are at the radiation-induced limit, while cases (c) and (d) are at the transport-induced limit. 

The unsteady simulations were carried out in order to determine the stability of the steady- 
state solutions. To do so, a steady-state temperature profile of the near- limit flames was perturbed 
and the ensuing transient response was noted. The perturbation either damped out and relaxed 
back to steady-state or triggered a growing oscillatory response resulting in extinction. 

Figure 3 shows the maximum temperature versus time for the oxidizer mixture given in case 
(a). For the 32.20% H 2 , the perturbed flame relaxes back to the steady-state solution. However, 
when the H 2 mole fraction is reduced to 32.19%, oscillatory instability develops and the flame 
eventually extinguishes after quite a number of cycles. The oscillation grows even faster for 
32.18%. Similar behavior was observed for case (b) except that the onset of oscillation was 
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28.09%. Figure 4 shows the maximum temperature versus time for cases (c) and (d). For case 
(c), pulsating extinction was not observed computationally; instead the perturbed flame at an 
ambient FF mole fraction of 8.024% would extinguish without any oscillations. For case (d), 
pulsating extinction was observed at H? = 7.727%, while FF = 7.728% is stable, although only 
one cycle of oscillation is observed before extinction. This also verifies that increase in Le in the 
oxidizer does make the flame pulsate before extinction. 

The experimental and computational results are compared in Table 1. It is seen that for 
oxidizer mixtures with CCF (cases (a) and (b)) the computational results over-predict the FF mole 
fraction at extinction by a substantial amount, while the opposite holds for cases (c) and (d). This 
is mainly due to the two heat transfer mechanisms that were neglected in the computations: (1) 
heat loss to the rod support and burner, which would promote extinction at higher FF mole 
fractions, (2) radiation reabsorption, which would strengthen the flame. Hence the CCF-diluted 
flames were experimentally observed to extinguish at lower fuel concentrations due to the 
dominance of reabsorption over heat loss to the support, while for cases (c) and (d), flames were 
observed to extinguish at higher fuel concentrations. Despite the disagreement on the hydrogen 
concentration at extinction, the experimentally observed frequencies at the onset of oscillation do 
agree quite well with the calculations, hence providing at least partial support to the notion that 
the same mechanism is responsible for both the computationally and experimentally observed 
oscillations. While pulsating extinction was observed computationally for case (d), pulsations 
were not observed experimentally. However, close inspection of the computationally observed 
pulsating extinction reveals that the flame extinguishes after only one cycle of oscillation, and the 
amplitude of oscillation is about 2K (compared to 10K for cases (a) and (b)), suggesting that 
these oscillations might prove difficult to be observed experimentally. 
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Maximum Temperature (K) Radiometer Reading (Volts) 


Table 1 Comparison of experimental and numerical results, where Leo is the oxidizer Lewis number and Le F the fuel 
Lewis number ( *No pulsation observed) 



Case (a) 

Case (b) 

Case (c) 

Case (d) 

Comp. 

Exp. 

Comp. 

Exp. 

Comp. 

Exp. 

Comp. 

Exp. 

Mixture Composition 

87.15% C0 2 / 
12.85% 0 2 

70.95% C0 2 / 16.20% He/ 
12.85% 0 2 

87.15% N2/ 
12.85% 0 2 

70.95% N 2 / 16.20% He/ 
12.85% 0 2 

Mass Flowrate (mg/s) 

59 

59 

39 

39 

%H 2 at Extinction 

32.19 

19.1 

28.09 

20.3 

8.024* 

15.9* 

7.727 

17.0* 

Lcf 

0.91 

1 

0.94 

0.99 

1.08 

1.02 

1.08 

1.014 

Le 0 

0.82 

0.82 

1.19 

1.19 

1.14 

1.14 

1.47 

1.47 

Frequency (Hz) 

2.3 

2.1 

2.3 

2.3 

n/a 

n/a 

0.1 

n/a 




0.020 

0.018 

0.016 
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0.012 

0.010 
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Figure 1 The radiometer voltage readings respect to time 
rear flame extinction for cases (a) and (b). For figure/case 
'a), steady state flame emission readings are also plotted at 
rydrogen molar fraction of 19.16%. 


Figure 2 The radiometer voltage readings respect to time near 
flame extinction for cases (c) and (d). 



Figure 3 Temporal response of the maximum temperature 
for case (a), with H 2 = 32.20, 32.19 and 32.18%. 



J i i U i I i i i £ I i i i i I i i i i I 

0 10 20 30 40 

Time (s) 


Figure 4 Temporal response of the maximum temperature 
for case (c), with H 2 = 8.025 and 8.024%, and case fd), 
with H 2 = 7.728 and 7.727%. 
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